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At a low resolution scale with Q 2 = fi 2 corresponding to the nucleon bound 
O state; deep inelastic unpolarized structure functions Fi(x,[i 2 ) and ^(x,^ 2 ) 



are derived with correct support using the symmetric part of the hadronic 
tensor under some simplifying assumptions in the Bjorken limit. For doing 



^ ■ this; the nucleon in its ground state has been represented by a suitably con- 



structed momentum wave packet of its valence quarks in their appropriate 
SU(6) spin flavor configuration with the momentum probability amplitude 
taken phenomenologically in reference to the independent quark model of 
scalar-vector harmonic potential. The valence quark distribution functions 
u v (x,/j?) and d v (x,n 2 ), extracted from the structure function F\{x,[i 2 ) in a 
parton model interpretation, satisfy normalization constraints as well as the 
momentum sum-rule requirements at a bound state scale of fi 2 = O.lGeV 2 . 
QCD evolution of these distribution functions taken as the inputs; yields at 
Ql = IhGeV 2 - Qq) and xd v {x,Q'o) in good qualitative agreement with 
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the experimental data. The gluon distribution G(x, Qq) and the sea-quark dis- 
tribution q s (x,Qo); which are dynamically generated using the leading order 
renormalization group equation; also match reasonably well with the available 

experimental data. 
12.39.pn,12.39.ki,13.60.Hb,13.90.+i 
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I. INTRODUCTION 



It is well known that low-energy description of hadron structure in terms of constituent 
quark models have been quite successful in explaining a large body of relevant experimental 
data. But at very high energies, quantum chromodynamics (QCD), the theory of strong 
interactions of quarks and gluons; sets a different framework of a more complex quark- 
parton picture of hadrons for understanding the deep-inelastic scattering(DIS) phenom- 
ena. In this picture the deep-inelastic lepton-nucleon scattering is described in terms of 
unpolarized structure functions Fi(x, Q 2 ) and F 2 (x, Q 2 ); which are expressed as the charge 
squared weighted combinations of quark-parton distribution functions f p (x,Q 2 ). These 
parton distribution functions f P (x,Q 2 ); interpreted as the probability of finding a parton 
p(quark or gluon) in the hadron with a fraction 'x' of the hadron momentum when probed 
with very high momentum transfer Q 2 ; play an important role in the standard model phe- 
nomenology providing a deeper understanding of the quark gluon structure of the hadron 
at very high energies. In this connection many experiments have been made to measure 
the deep-inelastic structure functions from which parton distributions inside the nucleon 
at very high energy have been extracted [1,2]. Although Q 2 -dependence of the parton dis- 
tribution functions(PDF) is successfully described by Dokshitzer-Gribov-Lipatov- Altarelli- 
Parisi(DGLAP) evolution equations [3] within perterbative QCD; absolute values of these 
observables are not provided theoretically by QCD to be compared with the experimental 
data. This is because; it requires some initial input distribution at lower resolution scale 
Q 2 = ji 2 which has not been possible from a first principle QCD-calculation due to the inad- 
equate understanding of the non-perturbative QCD in the confinement domain. Although 
lattice QCD as a favourite first principle technique has been pursued in this context [4]; 
it does involve inevitably increasing computational complexity in arriving at any desirable 
precision in its prediction. Therefore it had been a common practice to take the initial in- 
put distributions at a lower reference scale in suitable parametrized forms; which are fitted 
ultimately after the QCD evolution with the available experimental data. Alternatively, 
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there has been attempts to derive the distribution functions at the bound state scales of 
the nucleons described by the low energy QCD inspired phenomenological constituent quark 
models; which has been pursued over the years by many authors [5-13] with the purpose of 
establishing a much desired link between the low energy constituent quark picture and the 
high energy quark-parton picture of the hadron structure which may provide better under- 
standing of the parton distribution in nucleons inside the nucleus as well as of the parton 
contributions to the proton spin. 

The structure functions derivable from a constituent quark model corresponding to a low 
energy resolution scale Q 2 = /i 2 ~0(Aq CD )] is considered to represent the twist two non- 
singlet part of the physical structure function. Since at higher Q 2 -region; it is the twist two 
part of the physical structure function that dominates; QCD evolution of the model-derivable 
structure functions at Q 2 = ji 2 can provide results for comparison with the available data 
at higher-Q 2 . However the structure functions and the parton distributions derived at the 
bound state scale in constituent quark models usually encounter a pathological problem 
by not vanishing beyond x = 1 as required by energy-momentum conservation; which is 
commonly described as 'poor support'. Basing on the study of one-dimensional Bag model; 
Jaffe [14] had suggested a mapping of the distribution function so derived from the region 
< x < oo, to the kinematically allowed region < x < 1, which was applied to three 
dimensions as well, for removing the support problem. However this was just a prescription 
only. The problem has been addressed in the center of mass Bag model [7]; where an 
effective co-variant electromagnetic current of the nucleon is considered which satisfies the 
translational invariance and hence conserves the four momentum. Another approach of using 
the Pierels- Yoccoz projection was also suggested by Benesh and Miller [6]. Calculation based 
on Bethe-Salpter and light cone formalism [11] do avoid the support problem. Bickerstaff 
and Londergan [12] have tried with a different picture of the nucleon; where the confined 
constituent quarks are treated approximately as a system of infinite free fermion gas at finite 
temperature. Most of these early calculations with or without the support problem; yield 
more or less qualitatively reasonable results by way of fitting the experimental data with the 
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QCD evolved structure functions or the parton distributions realized from the model input 
expressions. 

In our earlier work,we also attempted to derive the structure functions of the nucleon 
at a low resolution scale in an alternative constituent quark model of relativistic indepen- 
dent quarks confined by an effective scalar-vector harmonic potential in a Dirac formalism; 
whose model parameters had been fixed earlier at the level of hadron spectroscopy and 
static hadron properties [15]. The predictive power of this model had also been successfully 
demonstrated in wide ranging low energy hadronic phenomena which include the weak and 
electromagnetic decays of light and heavy flavor mesons [16], elastic form factors and charge 
radii of nucleon [15], pion and kaon [17] and the electromagnetic polarizability of proton [18]. 
Extending this model to the study of deep-inelastic scattering of electrons off a nucleon; we 
had obtained quite encouraging predictions for the polarized structure functions g\ (x, Q 2 ) 
and g2(x,Q 2 ) [19] as well as the unpolarized structure functions F$(x, Q 2 ),F 2 n (a;, Q 2 ) with 
the resulting parton distributions [20] at a qualitative level. In these works we had taken the 
usual approximation that the nucleon at some static point of Q 2 = /i 2 ; consists only of the 
valence quarks with no gluons or sea-quarks as constituents. The model solutions for the 
bound valence quark eigen-modes provide the essential model input in expressing the elec- 
tromagnetic currents which ultimately define the relevant hadronic tensor for deep-inelastic 
process. Explicit functional forms of the polarized as well as unpolarized structure functions 
were then derived analytically from the antisymmetric and symmetric part respectively of 
the hadronic tensor in the Bjorken limit. However the structure functions so derived at 
the model scale expectedly encountered the support problem, although it was found to be 
minimal. Therefore in the present work; we would like to improve upon our earlier attempts 
by a somewhat different approach within the scope of the same model in order to realize 
correct support in the structure function from which the parton distributions in the nucleon 
can be extracted. 

For doing this; we describe the nucleon in its ground state by a suitably constructed 
momentum wave packet of its valence quarks in appropriate SU(6) spin-flavor configuration; 
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where each of these quarks is taken in its respective momentum states with a momentum 
probability amplitude derivable from its bound state energy eigen-mode obtained in the 
model. The wave packet includes explicitly a four delta function to ensure energy momen- 
tum conservation at the composite level. The quark field operators defining the electro- 
magnetic currents in the hadronic tensor are expressed as free field expansions. Then the 
unpolarized structure functions Ff(x,fi 2 ); derived from the symmetric part of the hadronic 
tensor with certain simplifying assumptions in the Bjorken limit, is found to be free from the 
support problem. It becomes also true for the valence quark distributions extracted from 
the structure function after appropriate comparison with its parton model interpretation, 
which furthermore satisfy the normalization requirements as well as momentum sum-rule 
constraints at the bound state scale. We therefore believe that these valence quark dis- 
tributions can provide adequate model based inputs for QCD-evolution to experimentally 
relevant higher Q 2 -region for a meaningful comparison with the experimental data. 

The paper is organized in the following manner. In sec-II; we discuss briefly the ba- 
sic formalism with necessary model inputs to describe the nucleon in its ground state as 
a wavepacket conserving energy-momentum from its constituent level of the three valence 
quarks taken in their respective definite momentum states with appropriate momentum 
probability amplitudes corresponding to their ground state eigen-modes. In sec-III we de- 
rive the unpolarized structure functions F^XjQ 2 ) and F 2 (x,Q 2 ) for the nucleon from the 
symmetric part of the hadronic tensor under certain siimplifying assumptions in the Bjorken 
limit. Sec- IV provides an appropriate parton-model interpretation of Fi(x, Q 2 ) leading to 
the extraction of the valence quarks distribution functions u v (x, Q 2 ) and d v (x,Q 2 ) at a 
model scale of low Q 2 = /i 2 . These valence distribution functions are found to satisfy the 
required normalization constraints. The bound state scale of Q 2 = fi 2 , which is not ex- 
plicitly manifested in the expressions for the distribution functions, is fixed on the basis of 
the renormalization group equations [13] by taking the experimental data of the momentum 
carried by the valence quarks at Qq = 15GeV 2 along with the same at Q 2 = fi 2 . The 
valence distribution functions u v (x, Q 2 ), d v (x,Q 2 ) are evolved to the higher reference scale 
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Ql = 15GeV 2 using the QCD non-singlet evolution equations; from which valence con- 
tributions to the structure functions such as [F^x, Qq)] v ,F2(x, Qq)] v and the combination 
\F% (x, Ql) — F£{x, Qq)] v = f[u v (x, Ql) — d v (x, Ql)} are evaluated for a comparison with the 
available experimental data. In sec-V; we attempt to obtain the gluon and the sea quark 
distributions G(x,Qq) and q s (x, Ql) respectively by dynamically generating them from the 
well known leading order renormalization group equation [21,22] with the valence distribu- 
tions as the inputs. Then we evaluate the momentum fraction carried by the quark sea, the 
gluons and the valence quarks at Ql = 15GeV 2 leading to the saturation of the momentum 
sum-rule. Finally to realize the complete structure functions F% ' n (x, Ql) and their difference 
[F%{x, Ql) — F^x, Ql)] taking into account appropriate sea contributions together with the 
corresponding valence parts; we consider some specific prescriptions for the flavor decom- 
position of the sea. The results are then compared with the available experimental data. At 
the end; sec- VI provides a brief summary and conclusion. 



II. MODEL FRAMEWORK 

In a parton model study of deep inelastic scattering (DIS) of electrons off the nucleon; 
which is pictured as three valence quarks embedded in sea of virtual quark antiquark pairs 
and gluons; the partons within the nucleon are treated as approximately free because of the 
asymptotic freedom property of QCD-interaction and light cone dominance of DIS. But from 
the point of view of a phenomenological quark model to start with; it may be quite justified 
to consider the nucleon as consisting only of three valence quarks, which eventhough might 
be dressed by the sea-quarks and the gluon; can be taken as the only resolvable individual 
units with no further discrenible internal structure at the hadronic scale of low Q 2 = fi 2 . 
The gluon and the sea-quark contents at Q 2 3> fi 2 ; can be realized through dynamic gen- 
eration via gluon bremsstrahlung and quark pair creation in the frame- work of QCD. The 
valence quarks constituting the nucleon at the model scale; being bound by the confining 
interaction within the hadronic volume; are not really free to be in any definite momentum 
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states. However in order to establish a link with the parton model picture of DIS; one can 
argue in principle that the bound valence quarks in a nucleon; during the virtual comp- 
ton scattering envisaged in the description of DIS; can be sensed by the interacting virtual 
photon in various momentum states with certain probabilities appropriate to their bound 
state energy-eigen modes. These momentum probability amplitudes can be realized from the 
fourier projections of their energy eigen modes. In that case the nucleon at the low resolution 
scale can be thought of as a bundle of free valence quarks in SU(6) spin flavor configurations 
with some appropriate momentum distribution satisfying in some heuristic manner the en- 
ergy momentum conservation. Then one can analyse the deep inelastic scattering in terms 
of free valence quarks interacting with the virtual photon at definite momentum states with 
specific momentum probabilities, which can enable one to establish a link between the low 
energy description of DIS with the parton model interpretation at high energy. 

In view of our above motivation; we prefer to represent the nucleon in its ground state 
with a definite momentum P and spin projection S, to a first approximation; by a normalized 
momentum wave packet of free valence quarks in the form; 

I p, s >= - r L= f {[^Lg^, k 2 , k 3 ) 

5\h + k 2 + k 3 - P) | T(k 1? k 2 , k 3 ; S) > (2.1) 

Here; | T(ki, k 2 , k 3 ; S) > provides the SU(6) spin flavor configuration of the valence quarks 
in definite momentum states expressed as; 

|T(k 1 ,k a ,k 8 ;5)> = £ ^^({A*} G S) 

l->-2,3 

x aJ 1 (k 1 ,Ai)a^(k 2 ,A 2 )aJ s (k 3 ,A 3 ) | > (2.2) 

we must mention here that ; Zq imm ({\} £ S) denotes the usual spin flavor co-efficients 
and a q (k, A), aj(k, A) are respectively the free quark annihilation and creation operators 
with definite momentum k and spin projection 'A'; which obey the usual anti commutation 
relations. Finally Gjv(ki, k 2 , k 3 ) represents the momentum profile function of the three 
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quarks which is subjected to the constraint of energy momentum conservation provided 
through the delta function in an adhoc manner. If we consider G q (k, A') as the momentum 
probability amplitude of the bound valence quark 'q' in its lowest energy eigenmode <3>+ A (r); 
to be found in a free state of definite momentum k and spin projection A'; then 

G s (p, A') = / dh *Wr) e ^-*pr) 



= G q {p)6rt (2.3) 

where; _E P = ^IpI 2 + m \ an d m <?(Pj A') is the usual free Dirac spinor. With reference to a 
specific phenomenological quark model such as the independent quark model with scalar 
vector harmonic potential [15], G q (k, A') can be worked out in the form [16] as; 

G q (k, A') = G q (k)5^ (2.4) 

when; 

G.(k) = (E * ± ± '"''I '^xp (-^g) (2.5) 



2 v / 2tt A 9 E p 
Here i2 9 is the ground state binding energy of the bound quark in the potential field V(r) = 
(1/2)(1 + 7°)(ar 2 + V ) with r 0g = (aA g )~^; and 

AT 2 = 8(£? g + m g ) 

Then the momentum profile function Gjv(ki, k 2 , k 3 ) of the three quarks in the nucleon can 
be expressed in the product form as ; 

G N (k lt k 2 , kg) = G(k 1 )G(k 2 )G(k 3 ) (2.7) 

It may however be noted that the momentum probability amplitude G ? (k) of individual 
quarks would be flavor independent in the non-strange sector; since the model adopted here 
assumes SU(2) flavor symmetry. Finally we have taken an overall normalization factor in 
Eq.(2.1) as [Af(P)}- 1 / 2 ] which can be determined considering the co-variant normalization 
condition, 
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<P,S\ P\ S' >= (2n) 3 2E N 5 3 (P - P')6 SS > (2.8) 

Using Eq.(2.1) in Eq.(2.8) and expressing the momentum probability distribution for quark 
qi as Pi(k) = |Gj(k)| 2 ; one can obtain, 

and; 

J n = / II ^P(k,)5 4 (A;i + h + k 3 - P) (2.10) 

The integral in Eq.(2.10) can be evaluated in a quark mass limit m q — >0 for the nucleon at 
rest. For doing this; we express the energy-delta function term appearing in the expression 

as; 



Jdkxl + |k 2 | + |(k x + k 2 )| — M) — y - ' * ' 2W 5(z - z) (2.11) 



(M — 


|ki 




|k 2 |) 






|k 2 | 





Here z = cos9 k2 ; which sets the limits of integrations for |k 2 | as ; 

M M 

(y - |ki|) < |k 2 | < — (2.12) 

Then defining; 

r— 

0y(|k|) = r d|k'| A (|k'|)p,(M-|k|-|k'|) 
J — l k l 

I N = / rf|k| Pl (|k|)p 23 (|k|) (2.13) 
Jo 

so that; the normalization constant for the nucleon state corresponding to its rest frame can 
be found as; 

^ P = °> = T§^ < 2 - 14 ' 

The integrals for ^(|k|) and I N can either be evaluated analytically or numerically. 
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III. STRUCTURE FUNCTIONS IN THE MODEL 



The hadronic tensor describing the deep-inelastic electron-nucleon scattering; which is 
expressed as the Fourier transform of single nucleon matrix element of the commutator of 
two electromagnetic currents in the form; 

Wvv = h I d ^ eiqi <p,s 1 J ^ (0)] 1 p,s> > (3,1) 

can be analysed in the present model frame-work to derive the nucleon structure functions. 
In Eq.(3.1) q is the virtual photon four- momentum and (P, S) are respectively the four- 
momentum and spin of the target nucleon, such that 

P^P^ = M 2 , S^S^ = -M 2 and P^S^ = (3.2) 

The coventional kinematic variables are usually defined as Q 2 = —q 2 > and x = Q 2 /2u; 
when v = P-q and 0<x<l. In the rest frame of the target nucleon; one takes P=(M, 0, 0, 0) 
and q=(v/M, 0, 0, ^v 2 /M 2 + Q 2 ). 

The hadronic tensor in Eq.(3.1) can be decomposed into a symmetric part Wffi and 
an antisymmetric part Wf$ respectively; when Wffl defines the spin averaged structure 
functions Fi(x,Q 2 ) and F 2 (x, Q 2 ) through a co-variant expansion in terms of the scalar 
functions Wi(x,Q 2 ) and W 2 {x,Q 2 ) as ; 

w<? = [ - g„ + "-f] HM*. Q 2 ) + [(P, - o^)(P, - *^)] ^jf 1 - (3.3) 

The deep-inelastic unpolarized structure functions F\(x, Q 2 ) and F 2 (x, Q 2 ) which become 
the scaling functions of the Bjorken variable x in the Bjorken limit (Q 2 — > oo, and v — > oo; 
with x fixed) are defined as F^x, Q^W^x, Q 2 ) and F 2 (x, Q 2 )=uW 2 {x, Q 2 )/M 2 . It is well 
known that while F 2 (x, Q 2 ) provides the contributions of the transverse virtual photons; 
a combination such as Wl(x,Q 2 ) = [F 2 (x,Q 2 )/2x — F±(x, Q 2 )} owes it to the longitudinal 
virtual photons. It can be shown that W L (x,Q 2 ) = ^^W^; so that with Wq as finite in 
the Bjorken limit; Wl — ► satisfies there-by the so called Callen Gross relation 
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F 2 (x,Q 2 ) = 2xF 1 (x,Q 2 ). (3.4) 

Now for a model derivation of the structure functions one can start with Eq. (3.1) with 
a static no gluon approximation for the target nucleon considered at rest with the nucleon 
state | P, S > represented as a momentum wave-packet of the constituent valence quarks as 
given in Eq.(2.9). However it is convenient to recast Eq.(3.1) into a more suitable form [5] 

as 

i w = 3»/^ <^te • ,o, /' ^3r / <^3^v ' , ^ r ' , 

x < P,S [^(r,t), J„(r',0)] \P,S>. (3.5) 

The electromagnetic current of the target nucleon is taken here in the form J M (£) = 

J2 e q4>q(£,) l^qiO > where e q is the electric charge of the valence quark of flavor q inside 
q 

the nucleon. The quark field operators Vv(0 i s expressed here appropriately as the free field 
expansion; 

MO = T^-rE / -&[a s (p,A)«(p,A) e -« + a;(p,A)T;(p,A) e i ^ (3.6) 
(2tt) 2 ±a J J2E P 



where the free Dirac spinors for the valence quarks taken in the zero mass limit as; 

u(p,\) = Je~ p ( )x\ 
\a-p/Ej 

v(p,\) = y/E p \ ^ jxx (3.7) 

Now expanding the curent commutator in Eq. (3.5) and taking the free quark propagator 
appearing in the expansion under impulse approximation written in the zero mass limit as; 

\imS D (x) = — ^ / d 4 kfbe{k )5{k 2 )e ±lkx , (3.8) 
m— >o [2it) i J 

where e(k ) = sign(k ); the symmetric part of the hadronic tensor Wffl can be obtained 

as; 

= [g,x9u, + g^gvx - g» v gx.]T x °. (3.9) 
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when 

T Act = [327r 4 5 3 (0)] _1 ^eJ / -^k x e(k )5(k 2 ) [ + °° dte^ +k ^ 

q J (27r) J -co 

x y rf 3 rrfVe- l(q+k) (r - r,) < A a >, (3.10) 

and 

< A CT >=< P,5 | [^(r,t) 7 ^,(r',0) -^(r',0) 7 ^(r,t)] | P, S > . (3.11) 

Since it is evident from Eq.(3.3) that Fi(x,Q 2 ) = Q 2 ) is the co-efficient of (—g^ u ) in 

the covariant expansion of W^ v \ Eq.(3.9) in the same token can yield; 

F l (x,Q 2 )=g Xa T x ° (3.12) 

Thus we find; 

F x {x,f) = [32vrV(0)]-7 -^e(k )5(k 2 ) J dte<<"+*<>>« 

x J d 3 rd 3 rV i(q+kHrV) < V >, (3.13) 

where; 

< r >=< P,S\J2 e 2 q [%(r, t)t%(T, 0) - #,(r', 0)^(r, *)] I p > ^ > (3-14) 

Now substituting the nucleon state | P, 5* > as in Eq.(2.1),(2.2), along with the free field 
expansions of the field operators as in Eq.(3.6),(3.7), we can realize after some necessary 
algebra 

F,(x, Q 2 ) = [f + (x, Q 2 ) - f_(x, Q 2 )} (3.15) 

where; 

f±M) = i6^(o) 7 ^ k 2jkn 

x / ^pdk.l) [l + ^ + <2E^)i] 5 3 (k + qikO^.o + k ±E kl ) 



/r/ 3 k r/ 3 k 
* - 3 p(|k 2 |)p(|k 3 |)^ 3 (k 1 + k 2 + k 3 )5(P fcl + E k2 + - M) (3.16) 



-^ikoiwi 
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It is to be noted here that, with the SU(2) flavor symmetry assumed in the present model; the 
spin flavor sum of the square of the quark charges of each flavor weighted by the respective 
probability |>2^| 2 corresponding to its SU(6) configuration denoted here as < e 2 >; gets 
decoupled from the rest of the integrals after simplification. Then one can independently 
evaluate < e 2 > for the spin up proton target as 1 and the same for the neutron target as 
2/3. 

In order to be able to perform the /^-integration first amongst the nested integrals, 
we first make a reasonable approximation to extract 5(k + q ±E kl ) from within the ki- 
integration as 5(ko+qo±E) with E kl —E corresponding to the peak position of the momentum 
distribution p(k\) under the expressions for f±(x,Q 2 ). It would then imply k = — (q ±E); 
which are always negative in the Bjorken limit for / + & /_ respectively. Then putting 
K = q + k so that |K| = K > K m = (|q| — |k|); where K m can reasonably be assumed to 
be much less than (q , |q| & |k|) in the Bjorken limit. Now doing the /co-integration; we can 
write for the proton; 

Q2) = vMw) I iipki(cos6k) f ^ 4( * - |k| ± B) 



X 



/(7 3 k (7 3 k 
— A— ^p(|k 2 |)p(|k 3 |)5 3 (k 1 + k 2 + k 3 )6(E kl + E k2 + E k3 - M) (3.17) 

The delta function 5(q — \\t\±E) in Eq.(3.17) sets the value of |k| = k as k± = q ±E 
and 5 3 (K±k!) sets the struck quark momentum as k x = q=K. This now leads to certain 
kinematic relations relevant in further simplifying the expression f±(x,Q 2 ) in the Bjorken 
limit, which are as follows; 

K m = K ± (x) = \(E T Mx)\, 

cos9 K ~ (MxtE)/K, 
cos9 K cos9 k ~ -(Mx=fE)/K, 

d(co86 k )k%~KdK. (3.18) 
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Now using these kinematic relations to-gether with the same procedure as described in 
Eq.(2.11) and (2.12) and finally substituting W(0) as in Eq.(2.13) to (2.14) after the neces- 
sary simplifications; we get 

M r°° dK 

fU^P 2 ) = JT /_ „ — 2 p{K)~p{K)\K - K ± (x)} (3-19) 
41 N Jk±(x) K z 

where we have used K±(x) = (E =F Mx). It may be noted here that p{K) represents the 
effects of the spectator quarks. 

Thus using Eq.(3.19) in (3.15); we can obtain the structure function F[(x,Q 2 ) for the 
proton at its bound state scale. Similar calculation can lead to F™(x, Q 2 ) for the neutron; 
which would be gFf (x, Q 2 ) in the present model with SU(2) flavor symmetry. Since as usual 
it can be shown here that Wq Q {x,Q 2 ) is finite [24] in the Bjorken limit which would lead 
to Wl — >0 satisfying the Callen-Gross relation from which ,n (x, Q 2 ) can also be realized 
using the expressions derived for Ff' n (x ) Q 2 ). 

IV. VALENCE QUARK DISTRIBUTION FUNCTIONS 

In a parton picture, if we define the quark parton distribution functions in the (u,d) 
flavor sector inside the nucleon in the usual manner as a combination of valence and sea 
components, such as, u{x, Q 2 ) = u v (x, Q 2 ) + u s (x, Q 2 ) and d(x, Q 2 ) = d v (x, Q 2 ) + d s (x, Q 2 ) 
with the corresponding antiparton distributions defined accordingly; then 

Ff(x, Q 2 ) = l/18[{4u(a;, Q 2 ) + d(x, Q 2 )} + {4u(x, Q 2 ) + d(x, Q 2 )}] 

F?(x, Q 2 ) = l/18[{4rf(x, Q 2 ) + u{x, Q 2 )} + {Ad(x, Q 2 ) + u(x, Q 2 )}] (4.1) 

Now comparing expressions in Eq. (4.1) with Eq. (3.15) and attributing as usual for such 
models ; the negative part of the distributions in Eq.(3.15) to the anti-partons in Eq.(4.1); 
effective parton distributions can be identified [5] as; 

u(x, Q 2 ) = 2d(x, Q 2 ) = Afl (x, Q 2 ) 

u(x, Q 2 ) = 2d(x, Q 2 ) = -4f p (x, Q 2 ) (4.2) 
15 



It is to be noted here that the negative antiparton distributions so obtained at the model 
scale calculation, can be treated only as a model artifact which infact is encountered in all 
such constituent quark models [5]. This spurious contribution needs to be appropriately 
eliminated in extracting the valence quark distribution correctly from the effective parton 
distributions in Eq.(4.2). Thus keeping in mind that u v (x, Q 2 ) = = d v (x,Q 2 ) as per our 
initial assumption and considering the spurious parton and anti-parton sea to be symmet- 
ric^. e u s (x) = u s (x) = u(x) and d s (x) = d s (x) = d(x) etc ); we get the appropriate valence 
distributions as 

u v (x, Q 2 ) = 2d v (x, Q 2 ) = A[f + (x, Q 2 ) + f p (x, Q 2 )}. (4.3) 

Thus the valence quark distribution functions u v (x,Q 2 ) and d v (x,Q 2 ) can be extracted 
at a model scale of low Q 2 = fi 2 in terms of analytically obtained expressions f±(x, Q 2 ) as 
functions of the Bjorken variable x; which can be evaluated by taking the model parameters 
(a, Vo) and other relevant model quantities such as (m q , E q , r 0q etc) described in sec-II as 
per their values found in its earlier applications in Ref [15,16] such as; 

(a,V ) = (0M71QQGeV 3 ,-0.1375GeV), 
(m q = m u = m d , E q , r 0q ) = (OMGeV, 0A5129GeV, 3.35227Ge\/" 1 ) (4.4) 

However; in view of such a current quark mass limit adopted in the model applications earlier; 
we believe in the justification of making all our calculations meant for the ultimate Bjorken 
limit with m q — >0 on the grounds of derivational simplicity. In that case the corresponding 
model quantities E q and r 0q , relevant for our calculations, are not much different from those 
given in Eq.(4.4); since their values now would be; 

(E q , r 0q ) = (O.U90GeV, 3.37489Ge\/" 1 ) (4.5) 

with the same potential parameters (a, V ) as in Eq.(4.4). 

we take here the actual physical mass of the proton M = 0.940Gel / and E~0.l8GeV 
corresponding approximately to the peak position of the momentum distribution p(|ki|). 
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The distributions xu v (x, Q 2 ) and xd v (x,Q 2 ) are evaluated numerically as functions of x 
which are presented in Fig.l and Fig. 2 respectively showing correct support. It is found that 
these distribution functions for the valence quarks satisfy the normalization requirement as; 




while the total momentum carried by the valence quarks at this low reference scale comes 
out as; 

f 1 xluJx, Q 2 ) + d v (x, Q 2 )}dx = 0.994 ~ 1 (4.7) 
Jo 

Thus with a close consistency in the requirements of normalization and momentum sat- 
uration at the model scale; it can be justfied to use these valence distributions as ap- 
propriate model scale inputs for QCD-evolution to higher Q 2 . Realizing the valence 
distributions at experimentally relevant higher Q 2 -region through QCD evolution; one 
can further evaluate the valence parts of the structure functions such as F%(x, Q 2 )] v = 
Q 2 ) and [F 2 n (x, Q 2 )] Q 2 ) as well as the valence part of the combination 

{F£(x,Q 2 ) -F?(x,Q 2 )] v = lxu v (x,Q 2 ). 

However the model scale of low Q 2 = ji 2 is neither explicit in the derived expressions for 
the structure functions nor in the valence distributions u v (x, Q 2 ) and d v (x, Q 2 ). Therefore 
we need to first fix the model scale Q 2 = /x 2 ; with the help of the renormalization group 
equation [13], as per which; 

f = A|c^ L ; (4.8) 
where; L = [V n=2 (Q 2 ) /\/™= 2 (/i 2 )] 1 /^s 2 /n(-^) and V n=2 (Q 2 ) = L 1 dxx\u v (x,Q 2 ) + 

QCD 

d v (x,Q 2 )} as the momentum carried by the valence quarks at Q 2 . Now taking the experi- 
mental reference scale Q 2 = 15GeV 2 for which V rn=2 (15G'e\/ 2 )~0.4 [2,8] and \/" =2 (/i 2 )~l as 
in Eq.(4.7) together with Aq C d = 0.232Gel / and a]f s 2 = 32/81 for 3-active flavors; one can 
obtain fi 2 = O.lGeV 2 . If one believes that the perturbation theory still makes sense down to 
this model scale for which the relevant perturbative expansion parameter a s (fi 2 ) /2n is less 
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than one (~ 0.358), one can evolve the valence distributions u v (x, /i 2 ) = 2d v (x, /i 2 ) to higher 
Ql, where experimental data are available. In fact one does not have much choice here, 
because taking any higher model scale on adhoc basis would require a non-zero initial input 
sea quark and gluon constituents for which one does not have any dynamical information 
at such scale and hence it would complicate the picture. Therefore when a s ([i 2 )/2ir is well 
within the limit to justify the applicability of perturbative QCD at the leading order and 
further since non-singlet evolution is believed to converge very fast [23], to remain stable 
even for small values of Q 2 /Aq CD ; one may think of a reliable interpolation between the 
low model scale of Q 2 = fi 2 < O.lGeV 2 and the experimentally relevant higher Q 2 » fi 2 , 
if one does not insist upon quantitative precision. With such justification and belief many 
authors in past have used the choice of low Q 2 = /i 2 (for example; /i 2 = 0.063GeV 2 [9], 
0.068GeV 2 ,0.09GeV 2 [10] and 0.06GeV 2 [12]) as their static point for evolution. Infact the 
choice of low Q 2 = fi 2 = O.lGeV 2 in such models is linked with the initial sea and gluon 
distributions being taken approximately zero at the model scale. Following such arguments; 
we choose to evolve the valence distributions by the standard convolution tecnique based on 
nonsinglet evolution equations in leading order [3,23] from the static point of fx 2 = O.lGeV 2 
to Ql = 15GeV 2 for a comparison with the experimental data. Our results for xu v (x, Qq) 
and xd v (x,Ql) at Ql = 15GeV 2 are provided in Fig. 1 and Fig. 2 respectively along with 
the experimental data, which on comparison shows satisfactory agreement over the entire 
range < x < 1. The valence components of the structure functions such as [F$(x, Ql)] v 
and [Frfix, Ql)] v together with the valence part of the combination [F$ (x, Qq) — F£(x, Qq)] v 
calculated at Qq = 15GeV 2 , are also compared with the respective experimental data in 
Fig. [3,4 &,5] respectively. We find that the agreement with the data in all these cases 
are reasonably better in the region x > 0.2. This is because in the small x region; the sea 
contributions to the structure functions not included in the calculation so far; are significant 
enough to generate the appreciable departures from the data as observed here. 

Therefore for a complete description of the nucleon structure functions and hence the 
parton distributions in the nucleon; the valence contributions discussed above need to be 
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supplemented by the expected gluon and sea-quark contributions at high energies. 



V. GLUON AND SEA QUARK DISTRIBUTIONS 

The gluon and the sea-quark distributions at high energy inside the nucleon can be 
generated purely radiatively with appropriate input of the valence distributions, using the 
well known leading order renormalization group (R.G)-equations [22,23]. Considering that 
at higher energyheavier flavors may be excited above each flavor threshold, we define the 
total sea quark distribution here upto three flavors as; 

q s (x, Q 2 ) = 2[u s (x, Q 2 ) + d s (x, Q 2 ) + s s (x, Q 2 )} (5.1) 

and the gluon distribution by G(X, Q 2 ). Their moments q™{Q 2 ) and G n (Q 2 ) respectively can 
be obtained in terms of the corresponding moment V n (Q 2 ) of the input valence distributions 
V(x, Q 2 ) = [u v (x, Q 2 ) + d v (x, Q 2 )} according to the RG-equations such as: 

G n {Q 2 ) = [ an{l ~ an) Lf s {L, a - - L a "}]V n (Q 2 ), (5.2) 

q n s {Q 2 ) = [Lf s {a n L~ a - + (1 - a n )L, al - }}V n (Q 2 ) , (5.3) 

where the n-th moments of the functions A(x, Q 2 )={G(x, Q 2 ), q s (x : Q 2 ), V(x, Q 2 )} are 
defined as 

A n (Q 2 )= ( l dxx n ~ 1 A(x,Q 2 ), (5.4) 
Jo 

and the RG-exponents such as {a n , (3 n , a^ s , a!^} in the conventional notations are derivable 
for the n-th moment as per Ref [23]. Finally L = = ^ff/A^j which can also be 

expressed here in terms of the momentum carried by the valence quarks at Ql on the basis 
of the momentum saturation by valence quarks at the model scale Q 2 = fj 2 as; 

f 1 dxxV(x,Q 2 )r 1/anNs2 ; (5.5) 
Jo J 
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With 

a NS = 32/81 for three active flavors considered here, the value of L Q comes out as 
L ~ 9. Then calculating the appropriate RG-exponents as per Ref [23] for n = 2,4,6,8 
(higher moments being significantly smaller are not considered here) and the corresponding 
moments V n (Ql = 15GeV 2 ) from the evolved valence distribution u v (x, Ql) = 2d v (x,Ql) 
at Ql = 15GeV 2 ; we evaluate the respective moments G n (Q 2 ) ) and ?"(<5o) from Eqns (5.2) 
and (5.3). Then gluon and sea-quark distributions can be extracted by a matrix inversion 
technique with the help of simple parametric expressions taken for xG(x, Ql) and xq s (x, Ql) 
as; 

xG(x, Ql) = [axx 2 + a 2 x + a 3 + a 4 /y/x ] (5.6) 

xq s (x, Ql) = [hx 2 + b 2 x + b 3 + b 4 /^/x~ ] (5.7) 

The moments calculated from these parametric expressions would now provide a set of 
simultaneous equations for each set of parameters {a^} and {hi} separately. Solving these 
equations by matrix inversion method we arrive at the values of these parameters as: 

{a h i = 1,2,3, 4}=(-0.8659, 2.0447, -2.1086, 0.9223) 

{bi,i = 1,2,3, 4}=(-0.2229, 0.5093, -0.5007, 0.2123) (5.8) 

Thus we generate somewhat reasonable functional forms for xq s (x,Ql) and xG(x, Ql) 
at Ql = 15GeV 2 which are provided in Fig. 6 and Fig. 7 respectively in comparison with 
the available experimental data. We find the qualitative agreement with the data quite 
encouraging with almost vanishing contributions in both cases beyond x > 0.5. 

We find next the momentum distributions for different constituent partons at 
Ql = 15GeV 2 by calculating the second moments of the distribution functions 
u v (x, Ql), d v (x, Ql), q s (x, Ql) and G(x, Ql) respectively so as to obtain them as; 

u v (Ql) = 0.279, (0.275 ± 0.011) 
d v (Q 2 ) = 0.140, (0.116 ± 0.017) 
q s (Ql) = 0.106, (0.074 ±0.011) 

G(Ql) = 0.475, (0.535) (5.9) 
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For a comparison; the experimental values are shown within the brackets against the 
calculated values. We find that the parton distributions realized at a qualitative level in the 
model at Ql = 15GeV 2 ; saturate the momentum sum-rule. Finally to evaluate the complete 
structure functions F2 P ' n \x, Ql) by supplementing the respective valence components with 
the necessary sea-contributions; we consider a flavor decomposition of the net sea-quark 
distribution q s (x,Ql). With an old option of a complete symmetric sea in SU(3)-flavor 
sector; 

[Fi(x,Q 2 )] sea = [F?(x,Ql)] sea = ^xq s (x,Q 2 ) (5.10) 

However it has been almost established experimentally that the nucleon quark sea is 
flovor asymmetric both in SU(2) as well as SU(3) sector. Experimental violation of Gottfried 
sum rule [25] and more recent and precise asymmetry measurements in the Drell-Yan process 
with nucleon targets [26] have shown a strong x-dependence of the ratio [d s (x, Q 2 )/u s (x, Q 2 )] 
with d s (x, Q 2 ) > u s (x, Q 2 ) for x < 0.2 and d s (x, Q 2 ) running closer to u s (x, Q 2 ) for x > 0.2; 
whereas around x = 0.18; d s (x, Q 2 ) — 2u s {x, Q 2 ). Neutrino charm production experiment by 
CCFR collaboration [26] also provides evidence in favor of the relative abundance of strange 
to non-strange sea quarks in the nucleon measured by a factor ^= ^jjfr^f ]> = 0.477±0.063. 
Therefore keeping these experimental facts in mind; we make a reasonable choice for the 
flavor structure of the sea quark distribution as defined in Eq. (5.1) by taking; 

d s {x,Q 2 Q ) = 2u s (x,Ql) 
s s (x, Ql) = l -[u s {x, Ql) + d s (x, Ql)] (5.11) 

Then we find the sea contributions to the structure functions F2 P ' n \x, Ql) as; 

m^Qi)] sea = l^s(x,Q 2 ) 

1 3 

[F2^Q 2 o)] sea = ^s(x,Q 2 ) (5.12) 

Then the complete structure functions F^x, Ql) and F^i^x, Ql) with the valence and 
quark sea components taken together are calculated and shown in Fig. 3 and Fig.4 We find 
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that the overall qualitative agreement is reasonable for the region x > 0.04. We have also 
shown in Fig.5; the structure function combination \F% (x, Ql) — F£{x, Ql)] by taking into 
account the asymmetric sea contribution as in Eq. (5.12), which provides a relatively better 
agreement with the available experimental data taken over a Q 2 -range [2] . 

VI. SUMMARY AND CONCLUSION 

Starting with a constituent quark model of relativistic independent quarks in an effec- 
tive scalar-vector harmonic potential and representing the nucleon as a suitably constructed 
wave-packet of free valence quarks only of appropriate momentum probability amlpitudes 
corresponding to their respective bound-state eigen-modes; we have been able to analyti- 
cally derive the deep-inelastic unpolarized structure function Ff (x, Q 2 ) at the model scale 
of low Q 2 = /i 2 = O.lGeV 2 with correct support. The valence quark distributions u v (x,Q 2 ) 
and d v (x,Q 2 ) have been appropriately extracted taking the parton model interpretation of 
Ff(x,Q 2 ). The valence distributions u v (x,Q 2 ),d v (x,Q 2 ) satisfy the normalization require- 
ment as well as the momentum sum-rule constraints; providing there-by suitable low energy 
model inputs for QCD-evolution to experimentally relevant Q 2 = Ql = !5GeV 2 .The valence 
distributions in the form xu v (x, Ql) and xd v (x, Ql); the valence components [F2 ,n (x, Ql)] v as 
well as the valence part of the combination [F${x, Ql) — F£(x, Ql)] are then realized through 
the QCD-evolution at Ql = 15GeV 2 ; which compare reasonably well with the experimental 
data in the expected range of the Bjorken variable. 

The gluon distribution G(x, Ql) and the total sea-quark distribution q s (x, Ql) are dynam- 
ically generated from the renormalization group equations taking the moments of the valence 
quark distributions at Ql = 15GeV 2 as inputs. The results for xG(x, Ql) and xq s (x, Ql) find 
good agreement with the experimental data. Calculation of the constituent parton momenta 
also yields the momentum percentage in the valence-quark sector as ~27.9% and ~14.0% 
for the 'u' and 'd' flavor quarks respectively; whereas in the sea-quark and gluon sector 
we find the same to be ^47.5% and ^10.6% respectively saturating there-by the expected 
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momentum sum-rule. Incorporating the sea-quark contributions to the valence part of the 
structure functions; the complete unpolarized structure functions F$ (x, Ql), Fg(x, Qq) and 
the combination [F^x, Ql) — F£{x, Qq)} are obtained in reasonable agreement with the data 
in the region x > 0.1. 

There are of course various finer features of the nucleon structure functions to-gether with 
their behaviour near the region x = 0; which would be beyond the limit of this simplistic 
approach in the model to address. Nevertheless,within its limitations, the model is found to 
provide a simple parameter free analysis of the deep-inelastic unpolarized structure functions 
of the nucleon leading to the realization of its constituent parton distributions at Ql = 
15GeV 2 with an over- all qualitative agreement with experimental data. 
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FIG. 3. The calculated [F%(x, Q 2 )] va i at Q 2 = /x 2 = O.lGeV 2 (dotted line) and its QCD evolved 
result at Qq = 15GeV 2 (dashed line). F% '(x, Q 2 )(valence+asymmetric sea: solid line) in comparison 
with experimental data taken from R.G.Roberts. et al in Ref.[l]. 
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FIG. 4. The QCD evolved result for [F£(x, Q 2 )] va i at Q 2 = 15GeV 2 (dashed line) and F 2 n (x, Q 2 ) 
(valence+asymmetric sea: solid line) in comparison with data from R.G.Roberts. et al in Ref [1]. 
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FIG. 5. The QCD evolved result for [F$(x, Ql)-F%(x, Qf)] (dot-dashedline-valence only; solid 
line- valence+asymmetric sea) at Qq = 15GeV 2 compared with the data, (data is over the Q 2 -range 
of the experiments as per Ref.[2]) 
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FIG. 7. The dynamically generated xG(x, Q 2 ) (solid line) at Qq = 15GeV 2 , compared with the 
data from T.Sloan et al in Ref [2]. 
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